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ABSTRACT 

The Kilo-Degree Survey (KiDS) is an optical wide-held survey designed to map the mat¬ 
ter distribution in the Universe using weak gravitational lensing. In this paper, we use these 
data to measure the density profiles and masses of a sample of ~ 1400 spectroscopically 
identified galaxy groups and clusters from the Galaxy And Mass Assembly (GAMA) sur¬ 
vey. We detect a highly significant signal (signal-to-noise-ratio ~ 120), allowing us to study 
the properties of dark matter haloes over one and a half order of magnitude in mass, from 
M ~ 10 13 — 10 14 ' 5 /i _1 Mq. We interpret the results for various subsamples of groups using 
a halo model framework which accounts for the mis-centring of the Brightest Cluster Galaxy 
(used as the tracer of the group centre) with respect to the centre of the group’s dark matter 
halo. We find that the density profiles of the haloes are well described by an NFW profile 
with concentrations that agree with predictions from numerical simulations. In addition, we 
constrain scaling relations between the mass and a number of observable group properties. We 
find that the mass scales with the total r-band luminosity as a power-law with slope 1.16±0.13 
(1-sigma) and with the group velocity dispersion as a power-law with slope 1.89 ± 0.27 (1- 
sigma). Finally, we demonstrate the potential of weak lensing studies of groups to discriminate 
between models of baryonic feedback at group scales by comparing our results with the pre¬ 
dictions from the Cosmo-OverWhelmingly Large Simulations (Cosmo-OWLS) project, ruling 
out models without AGN feedback. 

Key words: Cosmology: dark matter; Galaxies: haloes, large scale structure of the Universe; 
Physical data and processes: gravitational lensing; Methods: statistical 
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1 INTRODUCTION 

Galaxy groups are the most common structures in the Universe, 
thus representing the typical environment in which galaxies are 
found. In fact, most galaxies are either part of a group or have been 
part of a group at a certain point in time (Eke et al. 2004). How¬ 
ever, group properties are not as well studied compared to those 
of more massive clusters of galaxies, or individual galaxies. This 
is because groups are difficult to identify due to the small num¬ 
ber of (bright) members. Identifying groups requires a sufficiently 
deep 1 spectroscopic survey with good spatial coverage, that is near 
100% complete. Even if a sample of groups is constructed, the typ¬ 
ically small number of members per group prevents reliable direct 
dynamical mass estimates (Carlberg et al. 2001; Robotham et al. 
2011). It is possible to derive ensemble averaged properties (e.g.. 
More et al. 2009b), but the interpretation ultimately relies on either 
a careful comparison to numerical simulations or an assumption of 
an underlying analytical model (e.g., More et al. 2011) . 

For clusters of galaxies, the temperature and luminosity of the 
hot X-ray emitting intracluster medium can be used to estimate 
masses under the assumption of hydrostatic equilibrium. Simula¬ 
tions (e.g., Rasia et al. 2006; Nagai et al. 2007) and observations 
(e.g., Mahdavi et al. 2013) indicate that the hydrostatic masses are 
biased somewhat low, due to bulk motions and non-thermal pres¬ 
sure support, but correlate well with the mass. In principle, it is 
possible to apply this technique to galaxy groups; however, this is 
observationally expensive given their faintness in X-rays, and con¬ 
sequently samples are generally small (e.g., Sun et al. 2009; Eck- 
miller et al. 2011; Kettula et al. 2013; Finoguenov et al. 2015; Pear¬ 
son et al. 2015) and typically limited to the more massive systems. 

Furthermore, given their lower masses and the correspond¬ 
ing lower gravitational binding energy, baryonic processes, such as 
feedback from star formation and active galactic nuclei (AGN) are 
expected to affect groups more than clusters (e.g., McCarthy et al. 
2010; Le Brun et al. 2014). This may lead to increased biases in the 
hydrostatic mass estimates. The mass distribution in galaxy groups 
is also important for predictions of the observed matter power spec¬ 
trum, and recent studies have highlighted that baryonic processes 
can lead to significant biases in cosmological parameter constraints 
from cosmic shear studies if left unaccounted for (e.g., van Daalen 
et al. 2011; Semboloni et al. 2011, 2013). 

The group environment also plays an important role in de¬ 
termining the observed properties of galaxies. For example, there 
is increasing evidence that star formation quenching happens in 
galaxy groups (Robotham et al. 2013; Wetzel et al. 2014), due to 
ram pressure stripping, mergers, or AGN jets in the centre of the 
halo (Dubois et al. 2013). The properties of galaxies and groups 
of galaxies correlate with properties of their host dark matter halo 
(Vale & Ostriker 2004; Moster et al. 2010; Behroozi et al. 2010; 
Moster et al. 2013), and the details of those correlations depend 
on the baryonic processes taking place inside the haloes (Le Brun 
et al. 2014). Hence, characterisation of these correlations is crucial 
to understand the effects of environment on galaxy evolution. 

The study of galaxy groups is thus of great interest, but con¬ 
straining models of galaxy evolution using galaxy groups requires 
both reliable and complete group catalogues over a relatively large 
part of the sky and unbiased measurements of their dark matter 
halo properties. In the past decade, several large galaxy surveys 


1 Fainter than the characteristic galaxy luminosity L* where the power-law 
form of the luminosity function cuts off 


have become available, and significant effort has been made to reli¬ 
ably identify bound structures and study their properties (Eke et al. 
2004; Gerke et al. 2005; Berlind et al. 2006; Brough et al. 2006; 
Knobel et al. 2009). In this paper, we use the group catalogue pre¬ 
sented in Robotham et al. (2011) (hereafter R+ll) based on the 
three equatorial fields of the spectroscopic Galaxy And Mass As¬ 
sembly survey (herafter GAMA, Driver et al. 2011). For the reasons 
outlined above, determining group masses using “traditional” tech¬ 
niques is difficult. Fortunately, weak gravitational lensing provides 
a direct way to probe the mass distribution of galaxy groups (e.g., 
Hoekstra et al. 2001; Parker et al. 2005; Leauthaud et al. 2010). 
It uses the tiny coherent distortions in the shapes of background 
galaxies caused by the deflection of light rays from foreground 
objects, in our case galaxy groups (e.g., Bartelmann & Schneider 
2001). Those distortions are directly proportional to the tidal field 
of the gravitational potential of the foreground lenses, hence allow¬ 
ing us to infer the properties of their dark matter haloes without 
assumptions about their dynamical status. The typical distortion in 
the shape of a background object caused by foreground galaxies 
is much smaller than its intrinsic ellipticity, preventing a precise 
mass determination for individual groups. Instead, we can only in¬ 
fer the ensemble averaged properties by averaging the shapes of 
many background galaxies around many foreground lenses, under 
the assumption that galaxies are randomly oriented in the Universe. 

The measurement of the lensing signal involves accurate shape 
estimates, which in turn require deep, high quality imaging data. 
The shape measurements presented in this paper are obtained from 
the ongoing Kilo-Degree Survey (KiDS; de Jong et al. 2015). KiDS 
is an optical imaging survey with the OmegaCAM wide-field im¬ 
ager (Kuijken 2011) on the VLT Survey Telescope (Capaccioli & 
Schipani 2011 ; de Jong et al. 2013) that will eventually cover 1500 
square degrees of the sky in 4 bands ( ugri ). Crucially, the survey 
region of GAMA fully overlaps with KiDS. The depth of the KiDS 
data and its exquisite image quality are ideal to use weak gravi¬ 
tational lensing as a technique to measure halo properties of the 
GAMA groups, such as their masses. This is the main focus of this 
paper, one of a set of articles about the gravitational lensing analy¬ 
sis of the first and second KiDS data releases (de Jong et al. 2015). 
Companion papers will present a detailed analysis of the properties 
of galaxies as a function of environment (van Uitert et al. in prep), 
the properties of satellite galaxies in groups (Sifon et al. 2015), as 
well as a technical description of the lensing and photometric red- 
shift measurements (Kuijken et al. 2015, K+15 hereafter). 

In the last decade, weak gravitational lensing analyses of large 
optical surveys have become a standard tool to measure average 
properties of dark matter haloes (Brainerd et al. 1996; Fischer et al. 
2000; Hoekstra 2004; Sheldon et al. 2004; Parker et al. 2005; Hey- 
mans et al. 2006; Mandelbaum et al. 2006; Johnston et al. 2007; 
Sheldon et al. 2009; van Uitert et al. 2011; Leauthaud et al. 2012a; 
Choi et al. 2012; Velander et al. 2014; Coupon et al. 2015; Hud¬ 
son et al. 2015). However, the interpretation of the stacked lensing 
signal of haloes with different properties is not trivial. Haloes with 
different masses are stacked together, and a simple fit of the sig¬ 
nal using some function describing an average halo profile, like a 
Navarro-Frenk-White profile (Navarro et al. 1995, hereafter NFW) 
, can provide biased measurements. A natural framework to de¬ 
scribe the statistical weak lensing signal is the so-called halo model 
(Cooray & Sheth 2002; van den Bosch et al. 2013). It provides a 
statistical description of the way observable galaxy properties cor¬ 
relate with the mass of dark matter haloes taking into account the 
halo mass function, the halo abundance and their large scale bias. 

The outline of this paper is as follows. In Section 2, we sum- 


marise the basics of weak lensing theory. We describe the data used 
in this work in Section 3, and we summarise the halo model frame¬ 
work in Section 4. In Section 5, we present our lensing measure¬ 
ments of the GAMA galaxy groups, and in Section 6, we derive 
scaling relations between lensing masses and optical properties of 
the groups. We conclude in Section 7. 

The relevant cosmological parameters entering in the calcula¬ 
tion of distances and in the halo model are taken from the Planck 
best fit cosmology (Planck Collaboration et al. 2014): = 0.315, 

Qa = 0.685, a 8 = 0.829, n B = 0.9603 and Q h h 2 = 0.02205. 
Throughout the paper we use M 200 as a measure for the masses of 
the groups as defined by 200 times the mean density (and corre¬ 
sponding radius, noted as it^oo)- 


2 STATISTICAL WEAK GRAVITATIONAL LENSING 

Gravitational lensing refers to the deflection of light rays from dis¬ 
tant objects due to the presence of matter along the line-of-sight. 
Overdense regions imprint coherent tangential distortions (shear) 
in the shape of background objects (hereafter sources). Galaxies 
form and reside in dark matter haloes, and as such, they are bi¬ 
ased tracers of overdense regions in the Universe. For this reason, 
one expects to find non-vanishing shear profiles around galaxies, 
with the strength of this signal being stronger for groups of galax¬ 
ies as they inhabit more massive haloes. This effect is stronger in 
the proximity of the centre of the overdensity and becomes weaker 
at larger distances. 

Unfortunately, the coherent distortion induced by the host halo 
of a single galaxy (or group of galaxies) is too weak to be de¬ 
tected. We therefore rely on a statistical approach in which many 
galaxies or groups that share similar observational properties are 
stacked together. Average halo properties (e.g. masses, density pro¬ 
files) are then inferred from the resulting high signal-to noise shear 
measurements. This technique is commonly referred to as ‘galaxy- 
galaxy lensing’, and it has become a standard approach for measur¬ 
ing masses of galaxies in a statistical sense. 

Given its statistical nature, galaxy-galaxy lensing can be 
viewed as a measurement of the cross-correlation of some baryonic 
tracer 5 g and the matter density field <5 m : 

£ g m(r) = (<5 g (x)<5 m (x + r)) x , (1) 

where r is the three-dimensional comoving separation. The Equa¬ 
tion above can be related to the projected matter surface density 
around galaxies via the Abel integral: 

E (R) = Pm [ ‘ [1 + ^gm(V« 2 +n 2 )] an, (2) 

Jo 

where R is the co-moving projected separation from the galaxy, 7r s 
the position of the source galaxy, p m is the mean density of the 
Universe and II is the line-of-sight separation. 2 Being sensitive to 
the density contrast, the shear is actually a measure of the excess 
surface density (ESD hereafter): 

AE (R) = E(sJ R) - E (R), (3) 


2 Here and throughout the paper we assume spherical symmetry. This as¬ 
sumption is justified in the context of this work since we measure the lens¬ 
ing signal from a stack of many different haloes with different shapes, which 
washes out any potential halo triaxiality. 
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where E(^ R) just follows from E (R) via 

R) = E(#) R! dR'. (4) 

The ESD can finally be related to the tangential shear distortion 7 t 
of background objects, which is the main lensing observable: 

AE(7?.) = 7t(.R)S cr , (5) 


where 

* c 2 DM 

^ cr 47 tG D(zi)D(zu 2 S ) ’ 


(6) 


is a geometrical factor accounting for the lensing efficiency. In the 
previous equation, D(zi) is the angular diameter distance to the 
lens, D(zi, z s ) the angular diameter distance between the lens and 
the source and D(z 3 ) the angular diameter distance to the source. 

In the limit of a single galaxy embedded in a halo of mass M, 
one can see that Equation 1 further simplifies because £ gm (r) be¬ 
comes the normalised matter overdensity profile around the centre 
of the galaxy. The stacking procedure builds upon this limiting case 
by performing a weighted average of such profiles accounting for 
the contribution from different haloes. This is best formulated in the 
context of the halo model of structure formation (see e.g. Cooray 
& Sheth 2002, van den Bosch et al. 2013), and for this reason, we 
will embed the whole analysis in this framework (see Section 4). In 
Section 3.3, we describe how the ESD profile is measured. 


3 DATA 

The data used in this paper are obtained from two surveys: the Kilo- 
Degree Survey (KiDS) and the Galaxy And Mass Assembly survey 
(GAMA). KiDS is an ongoing ESO optical imaging survey with 
the OmegaCAM wide-field imager on the VLT Survey Telescope 
(de long et al. 2013). When completed, it will cover two patches of 
the sky in four bands ( u,g,r, i), one in the Northern galactic cap 
and one in the South, adding up to a total area of 1500 square de¬ 
grees overlapping with the 2 degree Field Galaxy Redshift survey 
(2dFGRS herafter, Colless et al. 2001). With rest-frame magnitude 
limits (5 <j in a 2” aperture) of 24.3, 25.1, 24.9, and 23.8 in the u, 
g, r, and i bands, respectively, and better than 0.8 arcsec seeing in 
the r-band, KiDS was designed to create a combined data set that 
included good weak lensing shape measurements and good pho¬ 
tometric redhifts. This enables a wide range of science including 
cosmic shear ‘tomography’, galaxy-galaxy lensing and other weak 
lensing studies. 

In this paper, we present initial weak lensing results based on 
observations of 100 KiDS tiles, which have been covered in all four 
optical bands and released to ESO as part of the first and second 
‘KiDS-DRl/2’ data releases to the ESO community, as described 
in de Jong et al. (2015). The effective area after removing masks 
and overlaps between tiles is 68.5 square degrees 5 . 

In the equatorial region, the KiDS footprint overlaps with the 
footprint of the GAMA spectroscopic survey (Baldry et al. 2010; 
Robotham et al. 2010; Driver et al. 2011; Liske et al. 2015), carried 
out using the AAOmega multi-object spectrograph on the Anglo- 
Australian Telescope (AAT). The GAMA survey is highly complete 


3 A further 48 tiles from the KiDS-DRl/2, mostly in KiDS-South, were 
not used in this analysis since they do not overlap with GAMA. 
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Figure 1. KiDS-ESO-DRl/2 coverage of the three equatorial GAMA fields (G09, G12, G15). Each grey box corresponds to a single KiDS tile of 1 square 
degree. The black circles represent groups with Nf 0 f ^ 5 in the G 3 Cv7 catalogue (R+l 1). The size of the dots is proportional to the group apparent richness. 
The filled red circles indicate the groups used in this analysis. These are all groups either inside a KiDS field or whose centre is separated less than 2 h -1 Mpc 
from the centre of the closest KiDS field. 


Table 1. Summary of the area overlap of KiDS-DRl/2 in the three GAMA 
fields and the number of groups with at least 5 members used in this anal¬ 
ysis. In parenthesis we quote the effective area, accounting for masks, used 
in this work. 


GAMA field 

KiDS-DRl/2 overlap (deg 2 ) 

Number of groups 

G09 

44.0 (28.5) 

596 

G12 

36.0 (25.0) 

509 

G15 

20.0 (15.0) 

308 


down to petrosian r-band magnitude 19.8 4 , and it covers ~ 180 
square degrees in the equatorial region, which allows for the iden¬ 
tification of a large number of galaxy groups. 

Figure 1 shows the KiDS-DRl/2 coverage of the G09, G12 
and G15 GAMA fields. We also show the spatial distribution of the 
galaxy groups in the three GAMA fields (open black circles) and 
the selection of groups entering in this analysis (red closed circles). 

Table 1 lists the overlap between KiDS-DRl/2 and GAMA 
and the total number of groups used in this analysis. Figure 2 
shows the redshift distribution of the GAMA groups used in this 
work and of the KiDS source galaxies, computed as a weighted 
sum of the posterior photometric redshift distribution as provided 
by BPZ (Benitez 2000). The weight comes from the lensfii code, 
which is used to measure the shape of the objects (Miller et al. 
2007) (see Sec. 3.2.1). The median redshift of the GAMA groups 
is z=0.2, while the weighted median redshift of KiDS is 0.53. The 
multiple peaks in the redshift distribution of the KiDS sources re¬ 
sult from degeneracies in the photometric redshift solution. This 
is dicussed further in K+15. The different redshift distributions of 
the two surveys are ideal for a weak lensing study of the GAMA 
groups using the KiDS galaxies as background sources. 

3.1 Lenses: GAMA Groups 

One of the main products of the GAMA survey is a group cat¬ 
alogue, G 3 C (R+ll), of which we use the internal version 7. It 
consists of 23,838 galaxy groups identified in the GAMA equato¬ 
rial regions (G09, G12, G15), with over 70,000 group members. 
It has been constructed employing spatial and spectroscopic red¬ 
shift information (Baldry et al. 2014) of all the galaxies targeted by 
GAMA in the three equatorial regions. The groups are found using 
a friends-of-friends algorithm, which links galaxies based on their 



Figure 2. Redshift distribution of the GAMA groups used in this analy¬ 
sis (red histogram) and the KiDS galaxies (blue lines). In the case of the 
GAMA groups, we use the spectroscopic redshift of the groups with at least 
5 members (R+ll), while for the KiDS galaxies the redshift distribution is 
computed as a weighted sum of the posterior photometric redshift distribu¬ 
tion as provided by BPZ (Benitez 2000). The weight comes from lensfd , 
used to measure the shape of the objects (Miller et al. 2007). The two verti¬ 
cal lines show the median of the redshift distribution of the GAMA groups 
and of the KiDS sources. The two peaks in the redshift distribution of the 
GAMA groups are physical (and not caused by incompleteness), due to the 
clustering of galaxies in the GAMA equatorial fields. 

projected and line-of-sight proximity. The choice of the linking 
length has been optimally calibrated using mock data (Robotham 
et al. 2011; Merson et al. 2013) based on the Millennium simula¬ 
tion 0 (Springel et al. 2005b) and a semi-analytical galaxy forma¬ 
tion model (Bower et al. 2006). Running the final group selection 
algorithm on the mock catalogues shows that groups with at least 5 
GAMA galaxies are less affected by interlopers and have sufficient 
members for a velocity dispersion estimate (R+ll). For this rea¬ 
son we use only such groups in our analysis. This choice leaves us 
with 1413 groups, in KiDS-DRl/2, 11% of the full GAMA group 
catalogue. 

Figure 3 shows the distribution of the total group r-band lu¬ 
minosity as a function of the redshift of the group, the group ap¬ 
parent richness, which is the number of members brighter that 


4 The petrosian apparent magnitudes are measured from SDSS-DR7 and 

they include extinction corrections (Schlegel maps) 5 (Q m , fb,, Q \, h, erg, n s )=(0.25, 0.045, 0.75, 0.73, 0.9,1.0) 
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Figure 3. Total group r-band luminosity as a function of the redshift of the 
group. The size of the points is proportional to the group apparent richness 
and the colour of the points indicates the group velocity dispersion corrected 
for velocity uncertainty. The shape of the distribution is typical of a flux 
limited survey. 


r = 19.8, and the group velocity dispersion corrected for veloc¬ 
ity uncertainty, for this subsample. These group r-band luminos¬ 
ity values are calculated by summing the r-band luminosity of all 
galaxies belonging to a group and targeted by GAMA and they 
also include an estimate of the contribution from faint galaxies be¬ 
low the GAMA flux limit, as discussed in R+ll. This correction 
is typically very small, a few percent at low redshift and a fac¬ 
tor of a few at z ~ 0.5 since most of the luminosity comes from 
galaxies around M* — 5 log h ~ —20.44 (Loveday et al. 2012, 
2015), and most of the groups are sampled well below M*. Note 
that all absolute magnitudes and luminosities used in the paper are 
k-corrected and evolution corrected at redshift z = 0 (R+ll). The 
global k-correction used by R+ll is compatible with the median k- 
correction of the full GAMA (McNaught-Roberts et al. 2014, Fig. 1 
in the paper). 

All the stellar masses used in this work are taken from Taylor 
et al. (2011), who fitted Bruzual & Chariot (2003) synthetic stellar 
spectra to the broadband SDSS photometry assuming a Chabrier 
(2003 ) IMF and a Calzetti et al. (2000) dust law. 


3.2 Sources: KiDS galaxies 

We measure the gravitational lensing effect induced by the GAMA 
groups using galaxy images from KiDS. We refer to K+15 for a 
detailed description of the pipelines used to measure shapes and 
photometric redshifts for those objects. We briefly summarise here 
the aspects of the data processing most relevant for this analysis. 


Kitching et al. 2008; Miller et al. 2013). lensfit takes full account 
of the point-spread function in the individual (dithered) exposures 
and prior knowledge of the ellipticity and size distributions of faint 
galaxies, returning an ellipticity estimate for each galaxy as well as 
an inverse variance weight that is related to the uncertainty of the 
measurement. 

The average number density of galaxies with lensfit weight w 
larger than 0, and satisfying the photometric redshift cuts described 
in the next section, is 8.88 per square arcmin, corresponding to an 
effective number density: 

2! 

n e s = ^2 Wi (7) 

i 

of 4.48 galaxies per square arcmin, where A is the survey area and 
<j% = 0.065 is the intrinsic ellipticity variance. This is a measure¬ 
ment of the statistical power of the weak lensing data (see Chang 
et al. (2013) and K+15 for more details). 

It is well known that shape measurements for galaxies with 
low signal-to-noise ratio and small sizes tend to be biased (e.g., 
Melchior & Viola 2012; Refregier et al. 2012; Miller et al. 2013; Vi¬ 
ola et al. 2014). This ‘noise-bias’ stems from the non-linear trans¬ 
formations of the image pixels involved in the derivation of galaxy 
image shapes. It has the form of a multiplicative bias, and a cal¬ 
ibration of the shape measurements is typically required in order 
to get an unbiased shear estimator. In this paper, we use the same 
calibration that was determined in Miller et al. (2013). This cali¬ 
bration depends on the signal-to-noise and the size of the objects 
and needs to be applied, in an average sense, to the recovered shear 
field. In addition to this multiplicative bias, shape measurements 
can also be affected by an additive bias caused by a non-perfect 
PSF deconvolution, centroid bias and pixel level detector effects. 
This bias can be empirically quantified and corrected for directly 
from the data, using the residual average ellipticity over the survey 
area. More detail on these ~ 10 per cent bias corrections can be 
found in K+15. 

The analysis presented in this paper has been applied to four 
different ellipticity catalogues. Three of these catalogues were gen¬ 
erated by rescaling all the ellipticity measurements by some factors 
unknown to the team and chosen by a collegue. Matthias Bartel- 
mann G , external to the collaboration. The amplitude of the rescal¬ 
ing has been chosen such that the cosmological parameters derived 
from a cosmic shear analysis using the four blind catalogues would 
not differ more than 10-cr, where sigma is the error from the Planck 
cosmological papers. We refer to this procedure as blinding, and 
we have used it to mitigate confirmation bias in our data analy¬ 
sis. The authors asked our external to unblind the true shear cat¬ 
alogues only just before paper submission. The authors were not 
allowed to change any of the results after the unblinding, without 
documenting those changes. Whilst the shear was blind, we did not 
blind measurements of group properties, such as their luminosity, 
or measurements of the source photometric redshifts. 


3.2.1 Shape measurements 

All of our lensing measurements are derived from the r-band expo¬ 
sures in KiDS. This is the band with the highest image quality of 
the survey, as the queue-scheduling at the telescope ensures that ob¬ 
servations in this filter are taken in the best seeing conditions. The 
images are processed with the THELI pipeline, which has been op¬ 
timized for lensing applications (Erben et al. 2013), and ellipticities 
for the galaxies are derived using the lensfit code (Miller et al. 2007; 


3.2.2 Photometric redshift measurements 

The observable lensing distortion depends on the distances to the 
lens and source (Equation 6). Redshifts to the lenses are known 
from the GAMA spectroscopy, but for the sources we need to re¬ 
sort to photometric redshifts derived from the KiDS-ESO-DRl/2 
ugri images in the ESO data release. Processing and calibration 

6 bartelmann@uni-heidelberg.de 
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of these images is done using the Astro-WISE environment (Mc¬ 
Farland et al. 2013), and flux and colour measurements use the 
‘Gaussian Aperture and Photometry’ (GAaP) technique designed to 
correct aperture photometry for seeing differences (Kuijken 2008). 
These colours form the basis of the photometric redshift estimates, 
obtained with BPZ (Benitez 2000; Hildebrandt et al. 2012). After 
extensive tests, we reject galaxies whose photometric redshift pos¬ 
terior distribution p(z) peaks outside the range [0.005,1.2] (K+15). 
In what follows the p(z) for each source is used in the calculation 
of distances, and in particular in the calculation of the critical sur¬ 
face density (see Equation 6). K+15 show that if the peak of each 
source’s p(z ) had been used as the estimate of the redshift, the av¬ 
erage value of S cr and hence the average ESD would have been 
underestimated by ~ 10 %. 


3.3 Measurement of the stacked excess surface density profile 

The shape measurement algorithm used in this work, lensfit, pro¬ 
vides measurements of the galaxy ellipticities (ei, e 2 ) with respect 
to an equatorial coordinate system. 

For each source-lens pair we compute the tangential et and 
cross component e x of the source’s ellipticity around the position 
of the lens, 

e t \ /— cos(2</>) -sin(2 <j>)\ fe i\ 

e x / V sin(2</>) -cos(2 <j>)) \e 2 ) ’ 

where is the position angle of the source with respect to the lens. 
The average of the tangential ellipticity of a large number of galax¬ 
ies in the same area of the sky is an unbiased estimate of the shear. 
On the other hand, the average of the cross ellipticity over many 
sources should average to zero. For this reason, the cross elliptic¬ 
ity is commonly used as an estimator of possible systematics in the 
measurements. Each lens-source pair is then assigned a weight 

wis = w s E“ 2 , (9) 


which is the product of the lensfit weight w s assigned to the given 
source ellipticity, and a geometric term E cr which downweights 
lens-source pairs that are close in redshift and therefore less sensi¬ 
tive to lensing. We compute the ‘effective critical surface density’ 
for each pair from the spectroscopic redshift of the lens Zi and the 
full posterior redshift distribution of the source, p(z s ): 


~ -1 _ 47tG 
- — 3 - 


f 


Di(zi)Dis(zi,z s 

D a (z s ) 


-p(z a )dz s . 


GO) 


Finally, following Equation 5, we compute the ESD in bins of pro¬ 
jected distance R to the lenses: 


AE (R) 


( J2ls WisEtS cr | 1 

V Eis wu j 1 + K(R) ’ 


(ID 


where the sum is over all source-lens pairs in the distance bin, and 


K(R) = 


E i s wi s m s 
E is ™ is 


( 12 ) 


is an average correction to the ESD profile that has to be applied 
to correct for the multiplicative noise bias m in the lensfit shear 
estimates. Typically, the value of the K(R) correction is around 
0.1, largely independent of the scale at which it is computed. 

Figure 4 shows the stacked ESD profile for all groups either 
inside a KiDS field or whose centre is separated by less than 2 
/i _1 Mpc from the centre of the closest KiDS field. It shows a 
highly significant detection of the lensing signal (signal-to-noise 




Figure 4. Top panel: ESD profile measured from a stack of all GAMA 
groups with at least 5 members (black points). Here, we choose the Bright¬ 
est Cluster Galaxy as the group centre. The open white circle with dashed 
error bars indicates a negative AE. The dotted red line and the dash-dotted 
blue line show the best fits to the data of NFW (Navarro et al. 1995) 
and singular isothermal sphere profiles, respectively. Neither of the single¬ 
parameter models provides a good fit to the data, highlighting that com¬ 
plex modelling of the signal is required. Bottom panel: ESD profile, mul¬ 
tiplied by R to enhance features at large radii, measured from the cross¬ 
component of the ellipticities for these same groups (blue points) and mea¬ 
sured around random points using the same redshift distribution of the 
groups (red points). We only use measurements at scales outside the dashed 
areas for the rest of the paper. 


ratio ~ 120). We note that the signal-to-noise is very poor at scales 
smaller than 20 /i _1 kpc. This is due to the fact that many objects 
close to the group centres are blended, and lensfit assigns them a 
vanishing weight. We exclude those scales from any further analy¬ 
sis presented in this paper. 

For reference, we also show the best fit singular isothermal 
sphere (SIS) and NFW models to the stacked ESD signal. In the 
case of the NFW model, the halo concentration is fixed using the 
Duffy et al. (2008) mass-concentration relation. Neither of the two 
single-parameter models provides a good fit to the data (Xred > 
2.5), highlighting how a more complex modelling of the signal is 
required (see Section 4). 

Figure 4 also includes two tests for residual systematic errors 
in the data: the cross-component of the signal and the signal mea¬ 
sured around random points in the KiDS tiles. On scales larger than 
2 ft -1 Mpc, small but significant deviations are evident. We believe 
that one possible origin of the non-vanishing signal around random 
points at these scales is due to the incomplete azimuthal average 
of galaxy ellipticities, but we cannot exclude some large scale sys¬ 
tematics in the shear data. The current patchy coverage of lensing 
data complicates a detailed analysis and here we simply note that 
the effect is small (less than 10 percent of the signal at 2 /i _1 Mpc) 
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and exclude data on scales larger than 2 /i - 1 Mpc. Future analyses 
based on more uniform coverage of the GAMA area from the KiDS 
survey will need to address these potential issues. 

To summarise, in the rest of the paper we will use only pro¬ 
jected distances in the range (0.02 — 2) fi^Mpc. Both the cross¬ 
component of the shear and the signal around random points are 
consistent with a null-detection over these scales. 


3.4 Statistical error estimate 


In a stacking analysis with many foreground lenses, the ellipticity 
of any source galaxy can contribute to the AS i estimate in mul¬ 
tiple radial bins i of different lenses. We summarize here how we 
compute the resulting covariances between the ESD estimates AEi 
from the data. 

We start from Equation 11, which gives the expression for 
AS,;. For simplicity, we drop in what follows the noise bias correc¬ 
tion factor 1 + K (R) as it can be considered to have been absorbed 
in the effective critical density S cr . 

We first rearrange the sum in Equation 11 to separate the con¬ 
tributions from each source s, by summing first over all lenses l 
that project within the radial bin i from source s; for each source s 
we denote this set of lenses as i s . We can then rewrite Equation 11 
as 


AEi = 


y ' W s (ei sCsi -f- £2sSsi) 

Y. Z s i 


where C, S and Z are sums over the lenses 


(13) 


C si = Y ~ E c t,Is C0S(2 (j)ls) , 
iei s 


(14) 


comparison to the case in which all groups are stacked together' 
and compute the signal in 10 logarithmically spaced radial bins be¬ 
tween 20 fi _1 kpc and 2 /i _ 1 Mpc. This leads to an ESD covariance 
matrix with 55 independent entries, which can be constrained by 
the 100 KiDS tiles used in this analysis. The corresponding ma¬ 
trix is shown in Figure 5 together with the correlation matrix ob¬ 
tained from Equation 17. The small but significant correlation be¬ 
tween the largest-radial bins is a consequence of the survey edges. 
We further show the diagonal errors obtained with the two meth¬ 
ods, labelled Analytical and Bootstrap. Based on the work by Nor- 
berg et al. (2009), we might expect that the bootstrapping technique 
leads to somewhat larger error bars, although on larger scales this 
trend may be counteracted to some degree by the limited indepen¬ 
dence of our bootstrap regions. Elowever, the conclusions of Nor- 
berg et al. (2009) are based on an analysis of galaxy clustering, and 
a quantitative translation of their results to our galaxy-galaxy lens- 
ing measurements is not easy and beyond the scope of this work. 
The difference between the error estimates using these two inde¬ 
pendent methods is at most 10% at scales larger than 300 fi _ 1 kpc. 

Based on the results of this test, we consider the covariance 
matrix estimated from Equation 17 to be a fair estimation of the 
true covariance in the data, and we use it throughout the paper. In 
our likelihood analyses of various models for the data (see next sec¬ 
tion), we account for the covariance between the radial bins as well 
as between the different lens samples used to compute the stacked 
signal. We note that future analyses with greater statistical power, 
for example those based on the full KiDS and GAMA overlap, and 
studies focusing on larger scales than those considered in this anal¬ 
ysis, will need to properly evaluate the full covariance matrix that 
incorporates the cosmic variance contribution that is negligible in 
this work. 


S si = Y sin(20 u ), 

iei s 


(15) 


and 


= E^ 


-2 

CT,ls ' 


(16) 


Since each 6k s is an independent estimate of the shear field, where 
k=l,2, the ESD covariance between radial bins i and j can then be 
easily written as: 


0 _ (J t w a (CsiCsj + SsiSsj) 

13 ( T,s w » z si){y, a w s z aj ) 


(17) 


where of = 0.078 is the ellipticity dispersion weighted with the 
lensfit weight, for one component of the ellipticity. We compute 
this number from the whole KiDS-ESO-DRl/2 area. 

Equation (17) can be generalised to also compute the covari¬ 
ance between the ESD estimates for two different lens samples m 
and n: 


COV„ 


£ s 


( Csi,mCsj,n T Ssi,mSsj : n ) 


(£ s 


0(£ s 


0 


(18) 


by restricting the sums for the C, S and Z terms to lenses in the 
relevant samples. 

We test the accuracy of the above calculation, which doesn't 
account for cosmic variance, against the covariance matrix obtained 
via a bootstrapping technique. Specifically, we bootstrap the signal 
measured in each of the 1-square degree KiDS tiles. We limit the 


4 HALO MODEL 

In this Section, we describe the halo model (e.g. Seljak 2000; 
Cooray & Sheth 2002), which we use to provide a physical interpre¬ 
tation of our data. We closely follow the methodology introduced 
in van den Bosch et al. (2013) and successfully applied to SDSS 
galaxy-galaxy lensing data in Cacciato et al. (2013). 

This model provides the ideal framework to describe the sta¬ 
tistical weak lensing signal around galaxy groups. It is based on 
two main assumptions: 

(i) a statistical description of dark matter halo properties (i.e. 
their average density profile, their abundance and their large scale 
bias); 

(ii) a statistical description of the way galaxies with different 
observable properties populate dark matter haloes. 

As weak gravitational lensing is sensitive to the mass dis¬ 
tribution projected along the line-of-sight, the quantity of inter¬ 
est is the ESD profile, defined in Equation 3, which is related to 
the galaxy-matter cross correlation via Equation 2. Under the as¬ 
sumption that each galaxy group resides in a dark matter halo, its 
average AE (R, z ) profile can be computed using a statistical de¬ 
scription of how galaxies are distributed over dark matter haloes 
of different mass and how these haloes cluster. Specifically, it is 
fairly straightforward to obtain the two-point correlation function, 

7 If the signal is split further into several bins according to some prop¬ 
erty of the group, we expect the relative contribution from cosmic variance 
compared to the contribution from shape noise to be even lower. 
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Figure 5. Left panel: ESD correlation matrix between different radial bins estimated from the data. This matrix accounts for shape noise and the effect of the 
mask and is computed as described in Section 3.4. Middle panel: ESD correlation matrix between different radial bins estimated using a bootstrap technique. 
It accounts for cosmic variance as well as shape noise. Right panel: Comparison of the square root of the diagonal elements of the two covariance matrices as 
a function of distance from the group centre (here the BCG). Note the lower noise in the left-hand panel and the small but significant correlation between the 
largest-radial bins, which is a consequence of the many survey edges. 


£ gm (r, z), by Fourier transforming the galaxy-dark matter power- 
spectrum, P gm (/c, z), i.e. 

£gm(r,z) = ^ JJ° P g ra(k,zf-^J^-k 2 Ak, (19) 

with k the wavenumber, and the subscript ‘g’ and ‘m’ standing for 
‘galaxy’ and ‘matter’. 

In what follows, we will use the fact that, in Fourier space, the 
matter density profile of a halo of mass M at a redshift z can be 
described as M Uh(k\M,z), where M = 47r(200p)f?2oo/3> and 
Uh(k\M) is the Fourier transform of the normalized dark matter 
density profile of a halo of mass M 8 We do not explicitly model the 
baryonic matter density profile (Fedeli 2014) because, on the scales 
of interest, its effect on the lensing signal can be approximated as 
that of a point mass (see Section 4.1). Because the lensing signal 
is measured by stacking galaxy groups with observable property 
O g r p, on scales smaller than the typical extent of a group, we have 
P g m(k, z) = P grpm (fc, z), where 

Pg r p m (fc, z) = J V(M\0 S rp) Hm(k, M, z) AM , (20) 

and 

Hm(k,M,z) = —u h (k\M,z), (21) 

Pm 

with p m the co-moving matter density of the Universe. Throughout 
the paper, the subscript ‘grp’ stands for ‘galaxy group’. 

The function P(M\O glp ) is the probability that a group with 
observable property O grp resides in a halo of mass M. It reflects 
the halo occupation statistics and it can be written as: 

P(M|C9 grp )dM = HgrpiM, z) rih(M, z) AM . (22) 


Here, we have used 


Pgrp (M, z) 


(N)o grp (M) 

n g rp(Og rp , Z) 


(23) 


where, (N)a eTp (AT) is the average number of groups with observ¬ 
able property O gTp that reside in a halo of mass M. 

Note that rih(M, z) is the halo mass function (i.e. the number 
density of haloes as a function of their mass) and we use the analyt¬ 
ical function suggested in Tinker et al. (2008) as a fit to a numerical 


N-body simulation. Furthermore, the comoving number density of 
groups, ft grp, with the given observable property is defined as 

n gTp (O gT p,z) = J(N)o srp (M)n h (M,z)AM . (24) 

Note that in the expressions above we have assumed that we 
can correctly identify the centre of the galaxy group halo (e.g., from 
the position of the galaxy identified as the central in the GAMA 
group catalogue). In Section 4.1, we generalize this expression to 
allow for possible mis-centring of the central galaxy. 

Galaxy groups are not isolated, and on scales larger than the 
typical extent of a group, one expects a non-vanishing contribution 
to the power spectrum due to the presence of other haloes surround¬ 
ing the group. This term is usually referred to as the two-halo term 
(as opposed to the one-halo term described in Equation 20). One 
thus has: 

Pgm {k) = P grp m (fc) + Pgrp m (&) • (25) 

These terms can be written in compact form as 

P g ^pm(k, z) = J ’H g rp{k,M,z)'H ln (k,M,z) n h {M,z) AM, 

(26) 


P g rpmi,k,z) — J AMiHgrplk, M\, z) rih(Mi,z) 

Jam, Pm(fc, M 2 , z) 7ih(M 2 , z) Q(k\Mi, M 2 , z) . (27) 

The quantity Q(k\M\,M,,z) describes the power spec¬ 
trum of haloes of mass Mi and M 2 . In its simplest imple¬ 
mentation 9 , used throughout this paper, Q(k\Mi, M 2 , z) = 
bh(Mi, z)bh(M,, z)P hn (k, z), where b\,(M,z) is the halo bias 
function and P lm (fc, z) is the linear matter-matter power spectrum. 
We note that, in the literature, there exist various fitting functions 
to describe the mass dependence of the halo bias (see for example 
Sheth et al. 2001 ; Sheth & Tormen 1999; Tinker et al. 2010). These 
functions may exhibit differences of up to ~ 10% (e.g. Murray 
et al. 2013). However, a few points are worth a comment. 

First, the use of the fitting function from Tinker et al. (2010) 
is motivated by the use of a halo mass function calibrated over the 


8 We use M 200 masses for the groups throughout this paper, i.e. as defined 9 See, for example, van den Bosch et al. (2013) for a more refined descrip- 

by 200 times the mean density (and corresponding radius, noted as /Too )- lion of this term. 
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same numerical simulation. Second, the halo bias function enters in 
the galaxy-matter power spectrum only through the two-halo term 
and as part of an integral. Thus, especially because we will fit the 
ESD profiles only up to R = 2/i -1 Mpc, the uncertainty related 
to the halo bias function is much smaller than the statistical error 
associated to the observed signal. 


4.1 Model specifics 

The halo occupation statistics of galaxy groups are defined via the 
function {N)o grp (M), the average number of groups (with a given 
observable property O grp , such as a luminosity bin) as a function 
of halo mass M. Since the occupation function of groups as a func¬ 
tion of halo mass, N grp (M), is either zero or unity, one has that 
(N)o srp (M) is by construction confined between zero and unity. 
We model (N)o grp (M) as a log-normal characterized by a mean, 
log[M/ (h~ 1 M @ )], and a scatter a IogSl : 


(N)o sip (M) oc 


V2tv 

Gy 


-exp 


log M 


(log M- log M) 2 

2(j 2 - 

log M 


• (28) 


We caution the reader against over-interpreting the physical 
meaning of this scatter; this number mainly serves the purpose of 
assigning a distribution of masses around a mean value. 

Ideally, for each stack of the group ESD (in bins of group 
luminosity or total stellar mass) we wish to determine both these 
parameters, but to keep the number of fitting parameters low we 
assume here that cr logJ g is constant from bin to bin, with a flat 
prior 0.05 ^ °io g M ^ 1.5. This prior does not have any sta¬ 
tistical effect on the results and it only serves the purpose of 
avoiding numerical inaccuracies. There is evidence for an increase 
in this parameter with central galaxy luminosity or stellar mass, 
(e.g. More et al. 2009b,a, 2011), but these increases are mild, and 
satellite kinematics (e.g. More et al. 2011) support the assump¬ 
tion that cr log ^[ is roughly constant on massive group scales (i.e. 
log[M/(/ i _1 Mq)] > 13.0). We have verified that our assumption 
has no impact on our results in terms of either accuracy or precision 
by allowing a logI ^ t0 be different in each observable bin. 

For each given bin in an observable group property, one can 
define an effective mean halo mass, (M), as 


(M)o sip = j V(M\O slp ) M AM 

J(N)o srp (M)n h (M,z)MdM 

n grp (Clgrp, z) 

where z is the mean redshift of the groups in the bin under consider¬ 
ation, and we have made use of Equation (22) and (24). The effec¬ 
tive mean halo mass, (M)o , is therefore obtained as a weighted 
average where the weight is the multiplication of the halo occupa¬ 
tion statistics and the halo mass function. 

The dark matter density profile of a halo of mass M, 
p m (r\M), is assumed to follow a NFW functional form: 


p m (r | AT) = 


5p 


(r/r s )(l +r/r s ) 2 


(30) 


where r s is the scale radius and 8 is a dimensionless amplitude 
which can be expressed in terms of the halo concentration parame¬ 
ter c m = R. 2 oo/r s as 


200 _ cf _ 

3 ln(l -I- Cm) - Cm/(1 + Cm) ’ 


(31) 


where the concentration parameter, c m , scales with halo mass. Dif¬ 
ferent studies in the literature have proposed somewhat different 
fitting functions (e.g. Bullock et al. 2001; Eke et al. 2001; Maccio 
et al. 2008; Duffy et al. 2008; Klypin et al. 2011; Prada et al. 2012; 
Dutton & Maccio 2014) to describe the relation c m (AT, z). Over¬ 
all, these studies are in broad agreement but unfortunately have not 
converged to a robust unique prediction. Given that those fitting 
functions have been calibrated using numerical simulations with 
very different configurations (most notably different mass resolu¬ 
tions and cosmologies), it remains unclear how to properly account 
for the above mentioned discrepancies. As these fitting functions 
all predict a weak mass dependence, we decide to adopt an effec¬ 
tive concentration-halo mass relation that has the mass and redshift 
dependence proposed in Duffy et al. (2008) but with a rescalable 
normalization: 

Cm(M200,z) = f c X c£ Uffy (M200,2) 

= / c x 10.14 ( afigftr) ~° 081 (1 + -)" 101 . (32) 

Note that at z = 0.25, one has Cm 11 ® 7 « 5 for halo masses with 
log [M/(h~ 1 Mq)\ « 14.3. We leave f c free to vary within a flat 
uninformative prior 0.2 ^ / c < 5. 

The innermost part of a halo is arguably the site where a ‘cen¬ 
tral’ galaxy resides. The baryons that constitute the galaxy may be 
distributed according to different profiles depending on the physi¬ 
cal state (for example, exponential discs for stars and /3-profiles for 
hot gas, see Fedeli 2014). The lensing signal due to these different 
configurations could in principle be modelled to a certain level of 
sophistication (see Kobayashi et al. 2015). However, at the smallest 
scales of interest here 10 , those distributions might as well be ac¬ 
counted for by simply assuming a point mass, Mp. In the interest 
of simplicity, we assume that the stellar mass of the brightest clus¬ 
ter galaxy (M+ CG (Taylor et al. 2011) is a reliable proxy for the 
amount of mass in the innermost part of the halo. Specifically, we 
assume that 

M P = ApM+ ca , (33) 

where Ap is a free parameter, within a flat prior between 0.5 and 5. 

The adopted definition of centre may well differ from the true 
minimum of the gravitational potential well. Such a mis-centring 
of the ‘central’ galaxy is in fact seen in galaxy groups (see e.g. 
Skibba & Maccio 2011 and references therein). George et al. (2012) 
offer further independent support of such a mis-centring, Ending 
that massive central galaxies trace the centre of mass to less than 
75 kpc//i. 

We model this mis-centring in a statistical manner (see also 
Oguri & Takada 2011, Miyatake et al. 2015, More et al. 2015a 
and references therein). Specifically, we assume that the degree of 
mis-centring of the groups in three dimensions, A(M, z), is pro¬ 
portional to the halo scale radius r s , a function of halo mass and 
redshift, and parametrize the probability that a ‘central’ galaxy is 
mis-centred as p 0 a. This gives 


T-Lgrp (k, M, z ) 


(AQ 0grp (M) 

ttgipG) 


(1— PoB+PoB 


xe {-0.5k 2 (r s 1Z o ff) 2 ]j 


(34) 


Setting either p 0 g or TZ 0 g to zero implies that there is de facto no 
offset. We treat the two as free parameters in Section 5. The pa¬ 
rameter p 0 ff, being a probability, is bound between zero and unity. 


10 We fit the data in the range 0.02 < R/(h 1 Mpc) < 2.0 
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We apply a flat uniform prior to lZ 0 f[ G [0,1.5]. We note that this 
prior is very conservative, as according to George et al. (2012) and 
Skibba & Maccio (2011) the mis-centring is expected to be smaller 
than the scale radius of a group, for which TZ 0 ff = 1. 

In summary, the model parameter vector, is defined as A = 
(logM^, <r log M, /c Ap,p oS ,U 0 g) where i = l...N hins . Through¬ 
out the paper, we bin group observable properties in 6 bins. This 
leads to a 11 parameter model. We use Bayesian inference tech¬ 
niques to determine the posterior probability distribution P(X\V) 
of the model parameters given the data, V. According to Bayes’ 
theorem, 


P{\[D) oc P(V |A) P( A) oc exp 



(35) 


where P{V\ A) is the likelihood of the data given the model param¬ 
eters,assumed to be gaussian, and P( A) is the prior probability of 
these parameters. Here, 


X 2 (^) — [AE fci y — AEfcj] T (C 1 )kk l ,ji' [AE*/ — AE y,j'] , 

(36) 


where AE^j is the j’th radial bin of the observed stacked ESD 
for the groups in bin k, and AEtj is the corresponding model 
prediction. C is the full covariance matrix for the measurements, 
computed as detailed in Section 3.4. 

We sample the posterior distribution of our model parameters 
given the data using a Markov Chain Monte-Carlo (MCMC). In 
particular, we use 11 a proposal distribution that is a multi-variate 
Gaussian whose covariance is computed via a Fisher analysis run 
during the burn-in phase of the chain, set to 5000 model evalua¬ 
tions. 


5 DENSITY PROFILE OF GALAXY GROUPS 

We measure the ESD signal around each GAMA group with at least 
5 members in 10 logarithmically-spaced radial bins in the range 
20 /i _1 kpc to 2 /i -1 Mpc. We first assign errors to those mea¬ 
surements by propagating the shape noise on the tangential shear 
measurement in each radial bin. We divide the groups into 6 bins 
according to a given observable property, such as their velocity 
dispersion, total r-band luminosity, apparent richness or r-band lu¬ 
minosity fraction of the BCG. Bin limits are chosen to make the 
signal-to-noise of the ESD roughly the same in each bin. Once the 
bin limits are defined, we compute the data covariance between ra¬ 
dial bins and between group bins as outlined in Section 3.4. We 
summarise the bin-limits, the number of groups in each bin, the 
mean redshift of the bin and the mean stellar mass of the BCG in 
Table 2 for the four observables considered in this work. 

The typical signal-to-noise ratio in each of the 6 luminosity 
bins is of order ~ 20-25. This is comparable to the signal-to-noise 
ratio reported by Sheldon et al. (2009) for a weak lensing analy¬ 
sis of ~ 130000 MaxBCG clusters using SDSS imaging, once we 
restrict the comparison to a similar luminosity range. 

We jointly fit the signal in the 6 bins using the halo model 
described in Section 4. Since GAMA is a flux limited survey, the 
redshift distributions of the groups in the six luminosity bins are 
different, as shown in Figure 6. When we fit the halo model to the 
data, we calculate the power spectra and mass function (Equations 
20-27) using the median of the redshift distribution in each bin. 

11 A python implementation of this sampling method is available via the 
MONTEPYTHON code thanks to the contribution by Surhud More. 


N 

(A 


il 

“ 2 L 0 )] 6(9.4,10.9] 

11. 

log[L grp /(tt- 2 L 0 )] 6(10.9,11.1] 

Jt 

loglLgip/tf 1 

,J 

- a L 0 )] 6(11.1,11.3] 

4 

log[L gIp /(ft" 2 L 0 )] 6(11.3,11.5] 

A. 


- a L 0 )] 6(11.5,11.7] 

1 

log[L grp /(A- 2 L 0 )] 6(11.7,12.7] 

JL 


0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

Z 


Figure 6. Redshift distributions of the GAMA groups used in this paper 
in the six r-band luminosity bins. The group luminosity increases from left 
to right and from top to bottom. The solid vertical black lines indicate the 
median of the distributions. 


For each observable property, we run 5 independent chains 
with different initial conditions. We evaluate the convergence of the 
MCMC by means of a Gelman Rubin test (Gelman & Rubin 1992), 
and we impose R < 1.03, where the R-metric is defined as the ratio 
of the variance of a parameter in the single chains to the variance of 
that parameter in an ” iiber-chain ”, obtained by combining 5 chains. 


5.1 Matter density profiles of group-scale haloes 

We first test whether the ESD measurements themselves support 
the halo model assumption that the group density profile can be de¬ 
scribed in terms of a mis-centred NFW profile with a contribution 
from a point-mass at small scales, and what constraints can be put 
on the model parameters. In the interest of being concise, we only 
present the results derived by binning the groups according to their 
total r-band luminosity (see Section 3), as statistically equivalent 
results are obtained when the groups are binned according to their 
velocity dispersion, apparent richness or r-band luminosity fraction 
of the BCG. The binning by other observables will become impor¬ 
tant in the study of scaling relation presented in Section 6. 

One needs to define the centre of the halo before stacking the 
ESD profiles of the groups. Following R+l 1, we have three choices 
for the group centre: the centre of light (Cen), the Brightest Clus¬ 
ter Galaxy (BCG) and the brightest galaxy left after iteratively re¬ 
moving the most distant galaxies from the group centre of light 
(IterCen). Throughout the paper, unless stated otherwise, we use 
the BCG as the definition of the centre, as it is a common choice in 
the literature. We investigate the effect of using the other two defi¬ 
nitions of the group centre in Section 5.1.4 and in Appendix A. 

Figure 7 shows the stacked ESD profiles (green points with er¬ 
ror bars) for the 6 bins in total r-band luminosity. Note that the error 
bars are the square root of the diagonal elements of the full covari¬ 
ance matrix, and we use dashed bars in the case of negative values 
of the ESD. The ESD profiles have high signal-to-noise through¬ 
out the range in total luminosity and in spatial scales. Red lines 
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Table 2. Summary of the bin limits used to compute the stacked ESD signal, the number of groups in each bin, the mean redshift of the groups in each bin and 


the mean stellar mass 

of the BCG. 




Observable 

Bin limits 

Number of lenses 

Mean redshift 

log«Mf CG [fe- 2 M 0 ])) 

l 0 g[Lgrp/(fo 2 L©)] 
o-/(s - 1 km) 

Nfof 

L BCG/ L grp 

(9.4,10.9,11.1,11.3,11.5,11.7,12.7) 
(0, 225, 325, 375, 466, 610,1500) 
(5,6,7,8,11,19,73) 

(1.0, 0.35, 0.25,0.18,0.13,0.08, 0) 

(540,259,178,233,142,66) 
(501,359,124,198,147,89) 
(481,261,170,239,181,86) 
(346,252, 296, 227,200,97) 

(0.13,0.20, 0.23, 0.26, 0.30, 0.35) 
(0.15,0.19, 0.21, 0.23, 0.26, 0.31) 
(0.21,0.21, 0.21, 0.19, 0.18, 0.16) 
(0.10,0.16, 0.20, 0.25, 0.29, 0.34) 

(11.00,11.23,11.29,11.37,11.47,11.70) 
(11.05,11.20,11.30,11.36,11.41,11.64) 
(11.17,11.23,11.29,11.29,11.35,11.45) 
(11.16,11.19,11.22,11.29,11.36,11.53) 



10' 1 10 ° 10 1 10 ° 10 1 10 ° 

Rf- 1 Mpc] 


Figure 7. Stacked ESD profile measured around the groups BCG of the 6 group luminosity bins as a function of distance from the group centre. The group 
l-band luminosity increases from left to right and from top to bottom. The stacking of the signal has been done using only groups with Nf 0 f ^ 5. The error 
bars on the stacked signal are computed as detailed in section 3.4 and we use dashed bars in the case of negative values of the ESD. The orange and yellow 
bands represent the 68 and 95 percentile of the model around the median, while the red line shows the best fit model. 


indicate the best-fit model, whereas orange and yellow bands indi¬ 
cate the 68 and 95% confidence interval. The model describes the 
data well with a reduced Xred = 1-10, 49 d.o.f, over the full scale 
range, for all the luminosity bins. This justifies our assumption that 
the ESD profile can be accurately modelled as a weighted stack of 
mis-centred NFW density profiles with a contribution from a point 
mass at the centre. 

The main results of this analysis can be summarised as follows 
(68% percent confidence limits quoted throughout): 

• For each r-band luminosity bin, we derive the probability that 
a group with that luminosity resides in a halo of mass M (see Equa¬ 
tion 22). We show the median of the probability distribution for the 
6 bins in Figure 8. We constrain the scatter in the mass at a fixed to¬ 
tal r-band luminosity to be c logI ^ = OYTlg'j’g. This sets the width 
of the log-normal distribution describing the halo occupation statis¬ 
tics. We remind the reader that cr log ^j is the width of the distribution 
in halo masses at given total luminosity of the groups and it is not 
the scatter in luminosity (or stellar mass) at a fixed halo mass that 
is often quoted in the literature and that one would expect to be 
considerably smaller (e.g. Yang et al. 2009; Cacciato et al. 2009; 
More et al. 2011; Feauthaud et al. 2012a). This hampers the pos¬ 
sibility of a one-to-one comparison with most studies in the litera¬ 
ture. However, we note that van den Bosch et al. (2007) and More 
et al. (2011) reported values of the scatter in halo mass at fixed 
luminosity that are as high as 0.7 at the bright end. Furthermore, 


More et al. (2015b) reported a value of 0.79lg;gg for the width 
of the low mass end distribution of the halo occupation statistics 
of massive CMASS galaxies. Given the non-negligible differences 
between the actual role of this parameter in all these studies, we 
find this level of agreement satisfactory. 

• For each luminosity bin, a mean halo mass is inferred with a 
typical uncertainty on the mean of ~ 0.12 dex. 

• The relative normalisation of the concentration-halo mass re¬ 
lation (see Equation 32) is constrained to be f c = 0.841q' 23, in 
agreement with the nominal value based on Duffy et al. (2008). 

• The probability of having an off-centred BCG is p 0 g < 0.97 
(2-sigma upper limit), whereas the average amount of mis-centring 
in terms of the halo scale radius, TZ 0 s, is unconstrained within the 
prior. 

• The amount of mass at the centre of the stack which con¬ 
tributes as a point mass to the ESD profiles is constrained to be 

M pm = A pm (Mf CG ) = 2.061^ {M? cg ). 


Figure 9 shows the posterior distributions of the halo model 
parameters and their mutual degeneracies. Table 3 and 4 list the me¬ 
dian values of the parameters of interest with errors derived from 
the 16th and 84th percentiles of the posterior distribution. We dis¬ 
cuss the constraints on the model parameters in further detail in the 
remainder of this section. 
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Figure 8. Probability that a group with a given r-band luminosity resides 
in a halo of mass M. The red lines show the median distribution, while the 
orange and the yellow contours show the 68 and 95 percentile around the 
median. 


5.1.1 Masses of dark matter haloes 


5.1.2 Concentration and mis-centring 

The shape of the ESD profile at scales smaller than ~ 200 /i _1 kpc 
contains information on the concentration of the halo and on the 
mis-centring of the BCG with respect to the true halo centre. How¬ 
ever, the relative normalisation of the concentration-halo mass re¬ 
lation, f c , and the two mis-centring parameters, p 0 ij and lZ 0 s are 
degenerate with each other. A small value of f c has a similar ef¬ 
fect on the stacked ESD as a large offset: both flatten the profile. 
To further illustrate this degeneracy, we show in Figure 11 the 2D 
posterior distribution of the average projected offset (p 0 ff x lZ a s ) 
and the normalisation of the concentration-halo mass relation. It is 
clear how a vanishing offset would require a low value of the con¬ 
centration. 

The derived constraints on the average projected BCG off¬ 
set are quite loose: p 0 u x lZ 0 s < 1.10r a (2-sigma). Hence one 
might argue in favour of a simpler model or a model with a less 
informative prior on IZofi. We address both aspects in the follow¬ 
ing ways. First, we run a version of the halo model on the same 6 
luminosity bins in which we assume no mis-centring (i.e. we as¬ 
sume that the BCG is always at the centre of the dark matter halo). 
We find a similar value of the reduced chi-squared (x? e d = 1-04, 
51 d.o.f.), comparable values for the 6 masses (always within one 
sigma) but tighter constraints for the relative normalisation of the 
concentration-halo mass relation, f c = 0.59lgJ 3 . This is perhaps 
not entirely surprising given that in this case f c is not degener¬ 
ate with any other model parameter. Second, we relax the prior for 
IZoft from 0 ^ lZ 0 fi ^ 1.5 to 0 ^ 1Z 0 a ^ 5. Also in this case, 
we find statistically equivalent halo masses and similar constraints 
on p 0 ff, lZ 0 ff, and f c as in the fiducial case. We summarise the 
results of these tests in Figure 12. We conclude that the fact that 
the reduced \ 2 values for the three model-configurations are very 
similar and always larger than unity suggests that the 11-parameter 
model is not too complex given the signal-to-noise of the data. Ig¬ 
noring the mis-centring in the model lowers the relative normalisa¬ 
tion of the concentration-halo mass relation to a 3-sigma deviation 
from the nominal value of Duffy et al. (2008). However, we caution 
the reader against over-interpreting this result as our test shows that 
this is probably driven by the very strong prior on the location of the 
BCGs rather than actually being a physical property of the stacked 
haloes. 


The dark matter halo masses of the galaxy groups that host the 
stacked galaxy groups analysed in this work span one and a half 
orders of magnitude with M £ [1O 13 ..1O 14 ' 5 ]/i _ 1 M0. Since our 
ESD profiles extend to large radii, our 2 /i _1 Mpc cut-off is larger 
than R 200 over this full mass range, these mass measurements are 
robust and direct as they do not require any extrapolation. The un¬ 
certainties on the masses are obtained after marginalising over the 
other model parameters. Typically these errors are 15% larger than 
what would be derived by fitting an NFW profile to the same data, 
ignoring the scatter in mass inside each luminosity bins. Note that 
a simple NFW fit to the data in the 6 luminosity bins, with fixed 
concentration (Duffy et al. 2008) would also lead to a bias in the 
inferred masses of approximately 25%. 

The inferred halo masses in each luminosity bin are slightly 
correlated due to the assumption that the scatter in halo mass is 
constant in different bins of total luminosity. We compute the cor¬ 
relation between the inferred halo masses from their posterior dis¬ 
tribution, and we show the results in Figure 10. Overall, the correla¬ 
tion is at most 20%, and this is accounted for when deriving scaling 
relations (see Section 6). 


Lower values of the normalisation of the concentration-halo 
mass relation from weak lensing analysis have been previously re¬ 
ported. For example, Mandelbaum et al. (2008) studied a sample 
of LRGs and MaxBCG clusters from SDSS and reported a 2-sigma 
deviation of the normalisation of the mass-concentration relation 
with respect to the simulation predictions. In this case, the lenses 
were assumed to be the true centre of the dark matter halo, and the 
analysis limited to scales larger than 0.5 /i _1 Mpc to limit the im¬ 
pact of mis-centring. From a weak lensing and clustering analysis 
of SDSS-III CMASS galaxies, Miyatake et al. (2015) also found 
a lower normalisation if mis-centring of the lenses is not included 
in the model but report agreement with the theoretical predictions 
once the mis-centring is included. A similar conclusion was de¬ 
rived by van Uitert et al. (2015) from a lensing analysis of LOWZ 
and CMASS LRGs from the Baryon Oscillation Spectroscopic Sur¬ 
vey (BOSS) SDSS-DR10 using imaging data from the second Red- 
sequence Cluster Survey (RCS2). In an analysis of the CFHT Stripe 
82 Survey for haloes of masses around 10 14 /i _1 Mq, Shan et al. 
(2015) also reported a nominal value of the normalisation of the 
concentration-halo mass relation lower than the Duffy et al. (2008) 
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prediction, but the discrepancy between observations and predic¬ 
tions from numerical simulations was not statistically significant. 

Possible explanations for a lower normalisation of the 
concentration-halo mass relation might include halo-triaxiality, 
which we do not account for in our model, substructures in¬ 
side the main halo (Giocoli et al. 2012), galaxy formation re¬ 
lated processes which can make halo density profiles shallower 
by expelling baryons into the outer region of the halo (Sales 
et al. 2010, van Daalen et al. 2011) and the assumed cosmologi¬ 
cal model. In fact, the value of the concentration at a given red- 
shift, as a measure of the formation time of haloes, depends on 
the background cosmology. To address this last point, we run the 
halo model assuming two alternative cosmologies: a slight devi¬ 
ation from the nominal Planck result (Q m , h,n s ,Qbh 2 ) = 
(0.302,0.818,0.68,0.9686,0.02197) (Spergel et al. 2015), and 
the best fit result of a clustering and lensing analysis on SDSS data 
{Qm,as,h,n s ,nbh 2 ) = (0.278,0.763,0.739,0.978,0.02279) 
(Cacciato et al. 2013), which we regard as an extreme change in 
light of the recent Planck results. We do not find any difference in 
the posterior distributions of any model parameters, in particular 
on f c . We hence conclude that our results are not affected by the 
assumed cosmology. 

5.1.3 Point mass: the innermost part of the halo 

Measurements of the ESD profile at scales smaller than ~ 50 
/i _1 kpc constrain the amount of mass at the centre of the halo. 
We model this as a simple point mass. The measured amplitude of 
the point mass is not degenerate with any of the other halo model 
parameters, demonstrating that, given the quality of the data, the 
details of the distribution of the baryons at the very centre of the 
haloes are not relevant to infer global properties of the dark matter 
halo, such as their masses or concentrations. 

5.1.4 Other definitions of the group centre 

Finally, we repeat the analysis using two alternative definitions of 
the group centre in the GAMA catalogue: the centre of light (Cen) 
and the brighter galaxy left after iteratively removing the most dis¬ 
tant galaxy from the group centre of light (IterCen). We present 
the results in Appendix A, Table 3 and 4. We do not find any signifi¬ 
cant difference in the ESD profile when using IterCen instead of 
the BCG. However, the profile is very different when we use Cen. 
In this case, we find tight constraints on the probability of the cen¬ 
tre of light of not being the centre of the dark matter halo with 
Pott ^ 0.67 at 2-sigma and we find that on average the amount 
of mis-centring of the centre of light with respect to the minimum 
of the halo potential well is lZ 0 s = I.OO^q ^). The constraints on 
the halo masses in the 6 luminosity bins, as well as the constraints 
on <T logJ ^, / c , and Ap, are however consistent within 1-sigma with 
those calculated using the BCG position. 

In summary, our results highlight the importance of a proper 
model for the mis-centring in the analysis of the ESD signal from 
groups or clusters of galaxies. Neglecting it could lead to biases 
in the derived parameters, particularly the normalisation of the 
concentration-mass relation. 


6 SCALING RELATIONS 

In the last Section of this paper, we investigate the correlations be¬ 
tween the halo masses derived using weak gravitational lensing and 



Bin 


Figure 10. Correlation matrix between the mean halo masses derived in 
the six r-band luminosity bins from the halo model fit. The reason for the 
correlation is the assumption of a constant scatter as a function of group 
luminosity in the halo occupation distribution. 



Figure 11. 2D posterior distribution of the average projected offset (p 0 fj X 
72 0 ff) and the normalisation of the concentration-halo mass relation f c . The 
contours indicate the 68%, 95% and 99% confidence region. 

optical properties of galaxy groups measured from SDSS images 
and the GAMA catalogue (R+ll). There are two main reasons to 
study these scaling relations: i) to understand which physical pro¬ 
cesses take place inside galaxy groups and their impact on galaxy 
formation; ii) to constrain a mean relation, as well as the scatter, be¬ 
tween some observable property of the groups and their halo mass 
for use in cosmological analyses that rely on the halo mass func¬ 
tion. 

6.1 The relation between halo mass and group r-band 
luminosity 

We first investigate the scaling relation between the total halo mass 
and the total r-band luminosity of the groups. As described in the 
previous section, we bin the groups according to their total r-band 
luminosity (see Table 2), fit a halo model to the stacked ESDs, and 
record the halo mass posteriors for each bin. We show the results. 
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Figure 9. Posterior distribution of the normalisation of the mass-concentration relation / c , of the mis-centring parameters p Q ff and 7£ 0 ff and of the amplitude 
of the point mass Ap. The contours indicate the 1, 2, 3 sigma confidence regions. The dashed vertical lines and the dotted vertical lines correspond respectively 
to the 1 and 2 sigma marginalised confidence limits. These are the constraints from a joint halo model fit of the ESD signal in the 6 luminosity bins using BCG 
as the group centre. The range in each panels reflect the priors used for the different parameters. 


Table 3. Constraints on the average halo mass in each r-band luminosity bin using the three definitions of halo centre. We quote here the median of the 
mass posterior distribution, marginalised over the other halo model parameters, and the errors are the 16th and 84th percentile of the distribution. All of the 
constraints derived using the three different proxies for the halo centre agree within 1-sigma. 


Centre 

loglM^Ah-iMo)] 

loglM'^/Oi-iMs)] 

log[M'| ) 0 /(h- 1 M o )] 

logKoo/( h- 1 M 0 )] 

log[M'| ) o /(h- 1 M 0 )] 

log[M' 6 0 ) 0 /(h- 1 M s )] 

BCG 

13.15±°;1| 

±o.oz _ 0 . 15 

13.83tS.i2 

i3.76tS:iS 

14 iq+°- 09 

14 ££+0.10 
i4.oo_ 0 10 

IterCen 

13.21±g'41 

13.45t° 0 iS 

13.76tSii 

13.77tSu 

14 16 +0 - 08 

14 kq+0.09 

J-4.04_q og 

Cen 

13.00tS.23 

13.64tS.i6 

13.92tS;iS 

13.85tS.ti 

14 1 q+0-09 

14.1o_q 10 

14 64+ 0 - 10 

14.O4_ 0 10 


P(A p ) 
















































































































KiDS+GAMA: properties of galaxy groups 15 


Table 4. Constraints on the halo model parameters using the three definitions of halo centre. For each of the parameters, we quote the median of the posterior 
distribution, marginalised over the other parameters, while the errors are the 16th and 84th percentile of the distribution. All the constraints derived using the 
three different proxies for the halo centre agree within 1-sigma. 
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Figure 12. Posterior distribution for the normalisation of the mass- 
concentration relation (Duffy et al. 2008) after marginalising over the other 
model parameters. We show here the effect of changing the prior in the mis- 
centring parameters: 7t 0 & = 0 (red line), TZ 0 g £ [0..1.5]r s (black line) 
and lZ 0 f{ E [0..5]r s (blue line). As a reference, the orange line shows the 
posterior distribution for / c in the case of a global stack of all groups. This 
has to be compared with the black line, where the constraints were derived 
from a joint fit of the stacked ESD in 6 luminosity bins. 


halo mass a function of group luminosity, in the left panel of Figure 
13. 

We fit a power-law relation between the halo mass and the 
total r-band luminosity of the group: 


M 200 

lO^/i-iMo 


(0.95 ±0.14) 


1O 11J5 /i- 2 L 0 


(1.16±0.13) 


(37) 


The linear regression is performed in the log-basis, since the errors 
on the masses are log-normal distributed, by minimizing the offset 
of the mass measurements from the power-law relation. We explic¬ 
itly account for the correlation between halo masses (see Section 
5). The red line in Figure 13 shows the best-fit relation. Our esti¬ 
mate of the 1-sigma dispersion around this relation is shown as the 
orange band and is derived from the joint posterior distributions for 
the halo masses from 5 independent MCMCs. We jointly extract 
10 5 random values of the masses in each of the 6 r-band luminos¬ 
ity bins (in order to preserve the correlation between the masses), 
and we fit a linear relation to each log-mass vector as a function 


of the logarithm of the r-band luminosity. Finally, we compute the 
16th and 84th percentiles of the best fit models in the different r- 
band luminosity bins. The average logarithmic scatter in halo mass 
at fixed r-band luminosity is ai og <m 20 o> = 0-1? 

In the left panel of Figure 13, we also compare our results to a 
previous weak lensing analysis of the same group catalogue (open 
black points) that used SDSS galaxies as background sources (Han 
et al. 2015). That analysis included all groups with Nf 0 f ^ 3 and 
fitted a single maximum likelihood mass to all the galaxies within 
a number of r-band luminosity bins. The agreement between the 
two analyses is remarkable given the different quality of data and 
the different techniques used to infer the halo masses. Neverthe¬ 
less, we stress that the current analysis based on the first KiDS data 
not only yields some of the tightest lensing constraints on group 
masses to date but also does this whilst marginalising over halo 
model parameters not considered in the previous work. 

Mock simulations suggest that the GAMA group catalogue 
is significantly contaminated by chance projections for groups 
with 2 and 3 members and marginally contaminated for groups 
with 4 members. Thus, while the only way to obtain constraints 
on low-luminosity systems (L grp < 10 w ' 5 Loh~ 2 ) is to include 
such sparse groups in the analysis, the impurity of the selection 
makes any results on the average mass of such groups unreliable 
and difficult to quantify (most likely underestimated). Our lowest- 
luminosity bin may suffer from a bias due to this same richness 
criterion if, as seems plausible, the poorer groups that are not in¬ 
cluded at a given luminosity have systematically lower masses. 

According to our current understanding of galaxy formation, 
one would expect the slope of the mass-luminosity relation to 
change towards the low-mass end, for haloes of about 10 12 — 
10 13 Mq/i _ 1 . This is mostly due to star formation being most effi¬ 
cient in haloes of ~ 10 12 /i _1 Mq (see for example Behroozi et al. 
2013 and references therein), implying the dominant feedback pro¬ 
cess is mass ejection from supernovae (see e.g. Dekel & Silk 1986). 
However, we are only able to probe the mass-luminosity relation 
for haloes more massive than about 10 13 Mgh -1 . In the regime 
modeled here, the relation is well fitted by a single power-law. 

The right panel of Figure 13 shows the relation between halo 
mass and total r-band luminosity in terms of the mass-to-light ratio. 
The mass-to-light ratio is relatively constant with total group r-band 
luminosity, with a slight increase of less than 0.1 dex from the low¬ 
est to the highest luminosity bin. The scatter around this ratio is 
as large as 0.2 dex. Ideally, one would like to compare this result 
with previous results from the literature. Unfortunately, different 
authors use different definitions of halo masses, group luminosities 
are often measured in different bands, and group selection functions 
might differ due to different survey depths or different algorithms 
used to identify groups. This might easily lead to different scaling 
relations, and we would like to highlight to the reader that a face- 
value comparison might be misleading. Despite these uncertainties. 
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Figure 13. Left panel: Halo mass as a function of the total group r-band luminosity. The solid black points show the halo masses derived in this work from 
a halo model lit to the stacked ESD profile of groups with at least 5 members brighter than the GAMA magnitude limit. The vertical error bars indicate the 
1-sigma uncertainty on the average halo mass after marginalising over the other halo model parameters, while the horizontal error bars indicate the 16th and 
84th percentile of the luminosity distribution in each bin. The red line shows the best fit power-law to the data points, while our estimate of the 1-sigma 
dispersion around this relation is shown as the orange area (see text). The open black circles show the halo masses derived from a lensing analysis of GAMA 
groups using SDSS galaxies as background sources (Han et al. 2015). Right panel: Derived mass-to-light ratio as a function of the group total luminosity from 
this work (black points), from the GAMA+SDSS analysis (open black circles), from the analysis of the CNOC2 group sample (Parker et al. 2005) (magenta 
diamonds) and from a lensing analysis of 130000 groups from the MaxBCG catalogue using SDSS imaging (Sheldon et al. 2009, (green crosses)). In blue we 
show the median relation derived using the 2PIGG catalogue (Eke et al. 2004). The red lines and the orange area correspond to those of the left panel. 


we qualitatively compare our results with previous measurements 
in what follows. 

One of the first analyses of a large sample of groups was 
based on the 2dFGRS, using a percolation technique to identify 
groups while also allowing dynamical mass measurements (Eke 
et al. 2004). The group luminosity was measured both in bj and 
r.F-band. We show this result as the blue line in the right panel 
of Figure 13. We find a qualitatively similar trend of the mass-to- 
light ratio as a function of the total group r-band luminosity for 
Lgrp > 10 11 L r@ h~ 2 . Flowever, our data do not support the steep 
increase of the mass-to-light ratio in the range 10 lo L r0 h~ 2 < 
L g ip < lO 11 L r0 h~ 2 reported by Eke et al. (2004). 

Han et al. (2015) carried out a detailed comparison between 
their results (which are in agreement with the one presented in this 
work) and the results from Eke et al. (2004), concluding that the 
steep increase in the mass-to-light ratio observed in the 2dFGRS 
sample could be mostly explained by the different depth between 
2dFGRS and GAMA (2 magnitudes deeper). We stress again here 
that our first data point might be affected by the apparent richness 
selection we applied on the group catalogue. If we exclude this data 
point, the agreement with Eke et al. (2004) is fairly reasonable. 

We also compare our results with a lensing analysis of 
MaxBCG clusters (Koester et al. 2007) using SDSS imaging (Shel¬ 
don et al. 2009). We show their result as the green points in Fig¬ 
ure 13. In this case the groups/clusters were binned according to 
their total luminosity and the masses were measured by first in¬ 
verting the ESD signal to 3D density and mass profiles and then 
by inferring the mass inside R 2 00 ■ Also in this case we find a rea¬ 


sonable agreement once we exclude our first data point, which, as 
discussed, might be affected by the apparent richness selection we 
applied to the group catalogue. 

Finally Parker et al. (2005) considered a sample of 116 groups 
from the CNOC2 survey (Yee et al. 1998). The halo masses were 
measured by fitting a SIS profile to the stacked ESD signal mea¬ 
sured using weak gravitational lensing. In this case, the luminosity 
was measured in i3-band. Given the small sample of groups, only 
two measurements were possible at quite low group luminosity. We 
show their results as the magenta points in Figure 13. Following lee 
et al. (2014), we applied a 0.8 multiplicative correction to the B- 
band mass-to-light ratio in order to have an estimate for the mass- 
to-light ratio in r-band. 

Only the mass-to-light ratio measurement in the high lumi¬ 
nosity bin of the CNOC2 analysis, which corresponds to our low 
luminosity bin, can be directly compared to our analysis, given the 
luminosity range we probe. We find a good agreement. 

6.2 The relation between halo mass and velocity dispersion 

Next, we focus on the scaling relation between the total halo mass 
and the group velocity dispersion. Again, we bin the groups in 6 
bins according to their velocity dispersion, with the boundaries cho¬ 
sen so that the signal-to-noise ratios of the stacked ESD profiles are 
equal (see Table 2). The halo masses in each bin are then found by 
a joint halo model fit to the ESD profile in each velocity dispersion 
bin. Figure 14 shows the corresponding results. The GAMA groups 
span an order of magnitude in velocity dispersion, but most of the 
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Figure 14. Halo mass as a function of the group velocity dispersion. The 
black points show the halo masses derived in this work from a halo model tit 
to the stacked ESD profile of groups with at least 5 members brighter than 
the GAMA magnitude limits. The red line shows the best fit power-law to 
the data points and the orange area indicates our estimate of the 1-sigma 
dispersion around this relation. The cyan points show the results from the 
CNOC2 survey (Carlberg et al. 2001), while the magenta points show the 
results from the HeCS sample of clusters (Rines et al. 2013). The grey band 
shows the mass-velocity dispersion relation obtained from measurements of 
satellite kinematics in SDSS (More et al. 2011). Finally, the blue line shows 
the relation calculated from the GAMA mocks using the same selection 
function applied to the data. 

constraining power for the scaling relation comes from groups with 
a ~ 500km s -1 . This is expected given that the cut imposed on 
group apparent richness excludes the low mass systems from this 
analysis, and that the survey volume is relatively small, and hence 
our sample does not contain many very massive galaxy clusters. As 
in the case of binning by luminosity, we believe that the apparent 
richness cut imposed on the GAMA group catalogue will have a 
non-neglible effect on the measurement of the average halo mass 
in the first velocity dispersion bin a < 200km s -1 . 

At low velocity dispersion, we compare our results with those 
from the CNOC2 survey (Carlberg et al. 2001), for which the mass 
measurements are derived from the dynamical properties of the 
groups. In Figure 14 we show the average CNOC2 mass measure¬ 
ments in 3 velocity dispersion bins; the error bars are the 1-sigma 
scatter between measurements in each bin. 

At high velocity dispersion, we compare our results to the 
analysis of the HeCS sample (Rines et al. 2013), where masses are 
measured using a redshift-space caustic technique. The mean red- 
shift of the HeCS clusters is similar to that of the GAMA groups. 
As for the CNOC2 sample, we binned the HeCS clusters accord¬ 
ing to their velocity dispersion, and we calculated the median mass 
and the 1-sigma dispersion in each bin. Both the CNOC2 and the 
HeCS sample agree well with the mass-velocity dispersion relation 
we derived using galaxy groups from GAMA. 

We fit a power-law between the halo mass and the group ve¬ 
locity dispersion (using the same procedure outlined in the previous 
section) and we constrain this relation to be: 

/ \ / N (l.S9±0.27) 

(lO^W) =^ 1 - 00±015 >(500s^k In ) 

( 38 ) 


We find that the average scatter in the halo mass-velocity dispersion 
relation is cr log<M200 > = 0.20. 

We do not see any indication of a change in the slope over 
almost two order of magnitude in mass, from massive clusters to 
small groups. However, the slope we find is significantly shallower 
than what would be expected from a virial scaling relation (M oc 
a 3 ) as is seen in dissipationless numerical simulation (Evrard et al. 
2008). A very similar result (M oc a 2 09± °' 34 ) was found by a 
previous weak lensing analysis of the same group catalogue using 
SDSS galaxies as background sources (Han et al. 2015). 

There are at least two possible explanations for this effect: 

• Hydrodynamical simulations have shown that galaxies trace 
shallower mass-velocity dispersion relations (slope lower than 3) 
than dark matter particles (Munari et al. 2013). This is due to dy¬ 
namical friction and tidal disruption, acting on substructures and 
galaxies, but not on dark matter particles. The typical effect mea¬ 
sured in simulations is of order 10%, which is too small to explain 
the value of the power-law slope we measure when comparing with 
the virial expectation. 

• The apparent richness cut we imposed to the group catalogue, 
the GAMA selection function and the limited cosmological vol¬ 
ume we probe might introduce selection biases on our mass mea¬ 
surements. In particular the apparent richness cut might introduce 
a positive bias for mass-measurements in the lowest velocity dis¬ 
persion bin, and the small volume used in this work might intro¬ 
duce negative biases in the highest velocity dispersion bins. The 
combination of these two effects would result in a shallower mass- 
velocity dispersion relation. 

To investigate the second hypothesis further we compare our in¬ 
ferred scaling relation with one measured from the dark matter only 
mock GAMA catalogue (Robotham et al. 2011; Merson et al. 2013) 
applying the same apparent richness cut. In the GAMA mocks the 
velocity dispersion is measured using the underlying/true dark mat¬ 
ter haloes while the stored mass of the haloes (DHalo mass) are 
computed as the sum of the masses of their component subhaloes 
(liang et al. 2014). For the purpose of the comparison we convert 
them into M 200 (McNaught-Roberts in prep.). We show the results 
as the blue line in Figure 14. We find a good agreement with the 
scaling relation measured from the data, supporting the hypothesis 
that the shallower scaling relation we measure is mostly caused by 
selection effects. However we cannot exclude at this stage that part 
of the reason for the shallower mass-velocity dispersion relation 
might be dynamical processes acting on the galaxies in the groups. 

A detailed investigation will be presented in a forthcoming 
paper (Robotham et al. in prep.) in the context of finding optimal 
dynamical mass estimates using weak lensing measurements of the 
group masses. 

Finally, we compare our results with measurements of the 
mass-velocity dispersion relation obtained from measurements of 
satellite kinematics in SDSS (More et al. 2011). In this case, we ex¬ 
trapolate the mass-velocity dispersion relation from measurements 
of the stellar mass - halo mass and stellar mass - velocity disper¬ 
sion relations which are provided in that paper. Note that these two 
relations have not been derived independently from each other. We 
find a good agreement with our results for a > 300kms 1 . For 
lower mass haloes, we have already discussed the potential selec¬ 
tion effect due to the apparent richness cut that affects our first data 
point. However, we also note that there is some tension between 
the CNOC2 results (Carlberg et al. 2001) and the SDSS satellite 
kinematics results. In general, velocity dispersion and mass mea¬ 
surements are more difficult for low mass groups than for massive 
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Figure 15. Group masses as a function of the r-band luminosity fraction of 
the BCG. The solid black points show the halo masses derived in this work 
from a halo model fit to the stacked ESD profile of groups with at least 5 
members brighter than the GAMA magnitude limit. The solid red and the 
dashed blue lines are predictions from the Cosmo-OWLS simulation at the 
median redshift of the GAMA groups for a run including AGN feedback 
and a reference run without AGN feedback (Le Bran et al. 2014). The lu¬ 
minosities measured in the simulation are (fc + e) corrected to redshift 0 
using the same functional form (Equation 39) applied to the data. The red 
area encompasses the 16th and 84th percentile of the mass distribution in 
each luminosity fraction bin for the AGN simulation. The shaded blue area 
indicates the range in Lbcg/ L grp in which there are no haloes in the REF 
simulation. 

systems because of the smaller number of members and more se¬ 
vere selection effects. 


6.3 The relation between halo mass and r-band luminosity 
fraction of the BCG 

Feedback from supemovae (Dekel & Silk 1986) and AGNs 
(Springel et al. 2005a) have been proposed in the past decade as 
a possible solution for reducing the star formation efficiency in 
hydrodynamical simulations (e.g. Sijacki et al. 2007, Fabjan et al. 
2010, McCarthy et al. 2010, Booth & Schaye 2013, Vogelsberger 
et al. 2014, Schaye et al. 2015 and references therein). It is impor¬ 
tant to test the hypothesis of feedback and to constrain its efficiency 
by comparing complementary predictions of hydrodynamical sim¬ 
ulations with observations. Motivated by the work of Le Brun et al. 
(2014), we focus here on the relation between the r-band luminos¬ 
ity fraction of the BCG, defined as Lbcg/L grp , and the group halo 
masses calculated in this work by binning the groups according 
to Lbcg /L grp (see Table 2). The r-band luminosity of the BCG 
is calculated from the tab petrosian magnitude from the GAMA 
catalogue. We apply a k-correction and evolution correction to the 
magnitude following R+11 : 

4 

(k + e){z) = Y, cii{z — 0.2)* — 1.75z , (39) 

i=0 

with a; = [0.2085,1.0226,0.5237, 3.5902, 2.3843]. We note that 
the original correction presented in Equation 8 in R+ll presents 
an error in the sign of the last term in the above equation. Figure 
15 shows the halo masses obtained for groups stacked according to 
Lbcg /L grp as a function of L B cG/L grp . 


There is a clear trend of group masses with the r-band lumi¬ 
nosity fraction of the BCG. This trend has been previously observed 
at group scales by Rasmussen & Ponman (2009) and at cluster 
scales by Lin & Mohr (2004). The explanation is that the growth of 
the BCG is modest compared with the growth of the entire group. 

Since the luminosity of the BCG is proportional to its stellar 
mass content and the group luminosity is an increasing function of 
the total halo mass, one can compare the results reported in this 
paper with studies of the stellar-to-halo mass ratio (SHMR) as a 
function of halo mass (e.g. George et al. 2011; Leauthaud et al. 
2012b; van der Burg et al. 2014; Coupon et al. 2015). There is 
a clear consensus on the decline of the SF1MR with halo mass, 
which is a different manifestation of the trend dispayed in Figure 
15 where we report the halo mass as a function of the r-band lumi¬ 
nosity fraction of the BCG. In particular, for central galaxies, it has 
been shown (Behroozi et al. 2013; Coupon et al. 2015) that halos 
of ~ 2 x 10 14 /i“ 1 Mq have a SF1MR about an order of magnitude 
lower than that of halos of ~ 10 13 /i _1 Mq, again in qualitative 
agreement with the result shown in Figure 15. The steep decline of 
the relation between the group mass and the r-band luminosity frac¬ 
tion is a consequence of star formation becoming less efficient in 
more massive halos. Several mechanisms, beyond AGN feedback, 
have been invoked to explain this phenomenon such as halo mass 
quenching (e.g Peng et al. 2010; Ilbert et al. 2013) or the presence 
of many satellite galaxies in massive halos which cut off the gas 
supply to the BCG (Aragon-Calvo et al. 2014). 

We focus here in particular on comparing our results with the 
(Cosmo-) Overwhelmingly Large Simulations (OWLS) (Schaye 
et al. 2010; Le Brun et al. 2014). 

Le Brun et al. (2014) present results from these simulations 
in terms of A'-band luminosity binned by halo mass. They report 
a very similar trend to the one we observe in our data. In par¬ 
ticular, they find a large difference in the luminosity fraction of 
the BCG when they compare simulations with and without AGN 
feedback. To compare our results with the Cosmo-OWLS sim¬ 
ulation, the r-band results were provided by the Cosmo-OWLS 
team using the same K-correction and evolution correction we 
applied to the data (Equation 39) for three redshifts snapshots 
2 = [0.125,0.25, 0.375], When comparing the simulations to the 
data, we use the results from the snapshots closer to the median 
redshift of the GAMA groups. We discarded from the simulation 
all haloes with mass lower than 10 13 /i _1 Mq, which roughly cor¬ 
responds to the minimum mass of groups with more than 5 mem¬ 
bers in the G 3 Cv7 catalogue (see section 6.4). In this way we try to 
mimic the selection we applied to the data. Finally, we bin the sim¬ 
ulation in the same way we bin the data, using the BCG luminosity 
fraction as a proxy for the group mass. 

Figure 15 shows the Cosmo-OWLS results for the run includ¬ 
ing AGN feedback (solid red line) and for a reference run (REF) 
without AGN (dashed blue lines). The red area encompasses the 
16th and 84th percentile of the mass distribution in each luminos¬ 
ity fraction bin. 

For LBCG/Lgr P < 0.2 the reference run does not contain any 
groups which on the contrary are clearly present in our group sam¬ 
ple. The reason for this is that the gas cooling in the REF simula¬ 
tion is too efficient, leading to BCGs which are always very lumi¬ 
nous in comparison to the total luminosity of the group. This evi¬ 
dence alone is sufficient to conclude that the data disfavour a model 
without AGN feedback. Note that this conclusion is independent 
of the group mass measurements. Our derived scaling relation be¬ 
tween the halo mass and the luminosity fraction of the BCG for 
Lbcg /L grp > 0.2 further supports the above conclusion, being in 
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reasonable agreement with the prediction from the simulation in¬ 
cluding AGN feedback. A detailed comparison of the trend in Fig¬ 
ure 15 with simulations would require replicating the GAMA group 
finder and selection function on the Cosmo-OWLS simulations and 
is beyond the scope of this paper. 

6.4 The relation between halo mass and group apparent 
richness 

Finally, we investigate the relation between the total halo mass and 
the apparent richness of the groups. The groups are binned accord¬ 
ing to their apparent richness (see Table 2), and the average halo 
mass for each bin is estimated by fitting a halo model to the stacked 
ESD profile. We show the result in Figure 16. 

We parametrise the halo mass-richness relation with a power- 
law, which is fit to the data with the same procedure outlined in the 
previous sections: 

= ^°' 43 ± °' 08 ^ \~W J ’ (40) 

and we constrain the average scatter in the halo mass-richness rela¬ 
tion to be rr log< M 2 oo> =0.20. 

As expected, richer groups are also more massive. We cau¬ 
tion the reader that this scaling relation is the one most affected 
by the GAMA selection function. In fact, unlike our treatment of 
the total group luminosity, we do not correct the apparent richness 
measurements to account for the faint galaxy members not targeted 
by GAMA. We compare our results with the GAMA mocks, which 
have the same selection function as the data, and we generally find 
good agreement. 

We also compare our results with a weak lensing analysis of 
130,000 groups and clusters of galaxies in the Sloan Digital Sky 
Survey (Johnston et al. 2007). The masses were derived from fit¬ 
ting an halo model to the stacked ESD profile in 12 richness bins. 
The richness was defined as the number of red sequence galaxies 
with luminosities larger than 0.4L* within a given projected radius, 
which is close to i? 200 - In spite of the different richness defintions 
we find a good agreement with our measurements, both for the am¬ 
plitude and the slope of the mass-richness relation. 


7 CONCLUSIONS 

In this paper, we present the first weak lensing analysis of the mass 
distribution in the GAMA groups using background sources from 
the overlapping KiDS survey. The effective overlapping area (ac¬ 
counting for masks) used in this work is 68.5 square degrees and 
corresponds to the first two data releases of ugri images of the 
KiDS data (de Jong et al. 2015 and K+15). 

Our main results are the following: 

(i) We measure the stacked excess surface density profile of the 
galaxy groups as a function of their total r-band luminosity, ve¬ 
locity dispersion, fraction of group light in the BCG and apparent 
richness. Splitting the data into six roughly equal signal-to-noise 
bins, we derive average halo masses per bin with a typical preci¬ 
sion of 0.12 dex. We provide a physical interpretation of the signal 
using the halo model. 

(ii) We show the importance of modelling the mis-centring of 
the BCG (used here as tracer of the group centre) with respect to the 
centre of the group’s dark matter halo in order to derive unbiased 
results, in particular on the halo mass-concentration relation. 



Figure 16. Group halo mass as a function of the richness (Nf 0 f) . We use 
here only groups with at least 5 members brighter than tab = 19.8. The 
richness of the groups is not corrected to account for the fainter galaxies not 
targeted by GAMA. The red line shows the best fit power-law relation to the 
data points. Our estimate of the 1-sigma dispersion around this relation is 
shown as the orange area. The blue line shows the mass-richness relation 
derived from the GAMA mocks using the same selection function applied 
to the data. The magenta points show the result of a weak lensing analysis 
of 130,000 groups and clusters of galaxies in the Sloan Digital Sky Survey 
(Johnston et al. 2007). 

(iii) Our results are consistent with the normalisation of the halo 
mass-concentration relation proposed by Duffy et al. (2008),when 
mis-centring is included in the model. 

(iv) We find no evidence of a significant baryonic component in 
the centre of the groups in excess of the stellar mass of the BCG. 
However, the uncertainty on this result is quite large due to the low 
signal-to-noise at small scales, which is in turn caused by the diffi¬ 
culties inherent in measuring reliable shapes for blended objects. 

(v) We obtain clear scaling relations between the halo mass and 
a number of observable properties of the groups: the group r-band 
luminosity, the velocity dispersion of the group, its apparent rich¬ 
ness and the ratio between the r-band luminosity of the BCG and 
the total r-band luminosity of the group. The typical scatter in halo 
mass at fixed observable property is r r iog(M 20 o) = 0-2- 

(vi) We show that our data have the statistical power to discrimi¬ 
nate between models with and without AGN feedback and possibly 
between different AGN feedback models. 

This analysis is part of the first set of weak lensing results 
using the KiDS data, based on data obtained during the first two 
years of operation. As the survey continues to cover more sky, both 
the statistical power and the fidelity of the measurements will grow, 
further refining these results as well as enabling other analyses of 
the distribution of dark matter in galaxies, groups and clusters. 
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APPENDIX A: ALTERNATIVE DEFINITIONS OF THE 
GROUP CENTRE 

We present here the measurements of the stacked ESD profile and 
the halo model constraints we obtain if we use a different definition 
of the group centre, compared to the BCG definitions used through¬ 
out the paper. 

In Figure Al we show the stacked ESD profile for the same 
6 luminosity bins used in Section 5 but now using the brightest 
galaxy left after iteratively removing the most distant galaxies from 
the group centre of light which is labelled as IterCen (left panel) 
and the group centre of light Cen (right panel) as the definition for 
the group centre. When IterCen is used, the stacked signal is sta¬ 
tistically indistinguishable from the case when BCG is used as the 
group centre. This is not surprising since the two centre definitions 
differ only for a few percent of the groups. 

When Cen is used as the group centre, the shape of the stacked 


ESD profile is very different. The turnover of the signal at scales 
around 100 /t _1 kpc is a clear indication of mis-centring between 
the chosen centre of the halo group and the true minimum of the 
halo potential well. R+l 1 report that Cen is not a good proxy for 
the halo centre, and hence, this result is not surprising. It is clear 
in this case that not including the mis-centring parameters in the 
model would lead to a very poor description of the data. 

We do not show the posterior distributions for the halo model 
parameters corresponding to the case of Cen. The degeneracies be¬ 
tween the parameters are the same as those found when BCG or 
IterCen are used as proxies for the halo centre. We can derive 
tight constraints on the probability of mis-centring p 0 ff ^ 0.67 
2-sigma, and we find that on average the amount of offset of the 
centre of light with respect to the minimum of the halo potential 
well is 77 0 ft = l.OOtolI- We summarise the results in Table 3 
and 4. 

The constraints we derive for the halo masses in the 6 luminos¬ 
ity bins and the constraints on cr log M, f c , Ap are consistent within 
1-sigma with the constraints derived using the other two definitions 
of the halo centre. 

These results highlight the importance of a proper model of 
the mis-centring in the analysis of the lensing signal from groups 
or clusters of galaxies. Neglecting mis-centring could lead to 
biases in the derived masses and in the other model parameters, 
particularly the halo concentration. 
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been performed considering only groups with Nf 0 f ^ 5. The error bars on the stacked signal are computed as detailed in section 3.4. The orange and yellow 
bands represent the 68 and 95 percentile of the model around the median and the red lines indicate the best fit model. 
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